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C/3 . The paper describes the transport of the elastically and diffractively 

^C ' scattered protons in the proton-proton interactions at the LHC for the 

• I high /3* runs. A parametrisation of the scattered proton transport through 

fi ■ the LHC magnetic lattice is presented. The accuracy of the unfolding of 

• ^ I the kinematic variables of the scattered protons is discussed. 
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1. Introduction 

> 

t:j- ' In higli energy proton-proton collisions at the LHC most attention is 

^ ■ usually paid to the hard processes. However, soft processes, such as elastic 

scattering or diffraction, contribute significantly to the total pp cross section. 
Studies of elastically scattered protons are important as this process can 
be used to precisely determine luminosity. Measurements of diffracively 
scattered protons can contribute to a better understanding of the still not 
well known soft QCD. 

Scattering angles of protons originating from elastic and diffractive in- 
teractions are very small, of the order of microradians. In order to reach 
this angular region in a collider environment, dedicated detectors must be 
C^ ' installed far away (dozens of metres) from the Interaction Point (IP) and 

very close to the beam (actually the detectors need to be placed inside the 
beam pipe). It is important to point out that typically there are several ac- 
celerator magnets between the IP and such detectors. Therefore, the proton 
trajectory depends not only on the scattering angle but also on the proton 
energy. It will be shown that from the measurement of the position and 
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direction of the proton trajectory one can obtain the information about all 
components of its momentum after the interaction. 

At the LHC there are two sets of detectors that are foreseen to measure 
elastically and diffractively scattered protons - the TOTEM experiment 
placed around the CMS Interaction Point and the ALFA detectors of the 
ATLAS experiment. This analysis focuses on the ALFA case. 

The ALFA experiment aims to determine the absolute luminosity of the 
LHC at the ATLAS from the measured rate of elastic scattering events in the 
Coulomb- nuclear amplitude interference region [Ij. It is worth mentioning 
that the ALFA detectors offer a possibility to study other processes, eg. 
single diffraction or even exclusive production [2j. For such measurements 
a special tune of the LHC machine is required. Such dedicated optics must 
deliver: 

(a) a very large value of the betatron function at the IP (/3*), 

(b) 90° phase advance of the betatron function between the IP and the 
detector locations in at least one transverse direction, 

(c) a small emittance (e*) of the beams. 

In fact, such optics provides parallel-to-point focusing in the (y, z) plane. A 
solution fulfilling the above requirements is called the high /3* optics. 

In the first phase (hereafter called the early high f3* ) the LHC is foreseen 
to run with the beam of 3.5 TeV energy, /3* = 90 m and e* = 2.5 /xm-rad. 
In the second phase (the nominal high j3*) the beam energy will be equal 
to 7 TeV, (3* to 2625 m and e* to 1 //m-rad. One should note that the 
nominal LHC (3* value for high luminosity runs is 0.55 m. The machine 
special settings for the high j3* optics are described in [3lll]. 

2. Transport Simulation 

To compute the particle trajectory in a magnetic structure of an acceler- 
ator one of several dedicated transport programs can be used. In this paper 
MAD-X [5], a program used to design and simulate particle beam behaviour 
within an accelerator, was employed. The program allows to perform the 
calculations using the thick lens approximation of the magnets - the Poly- 
morphic Tracking Code (PTC) module [6j. This module takes into account 
not only the magnetic structure and the geometry of the beam chamber, 
but also the fringe fields and edge effects. It is important to point out that, 
contrary to the studies of the beam, the thin lens approximation is valid 
for protons of the beam but not for the ones scattered in interactions. This 
is because the latter are more deflected in the magnetic field, hence their 
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distance from the magnets centres can be large and additional effects can 
play an important role. 

The LHC magnetic lattice in vicinity of the ATLAS IP is presented in 
Fig. [TJ Quadrupole magnets are labelled with letter Q and dipole magnets 
with letter D. In the following, a reference frame with the x~axis pointing 
towards the accelerator centre, the y-axis pointing upwards and the z-axis 
along one of the beams is used. All presented calculations were performed 
for the beam 1 that performs the clockwise motion. However, the results are 
qualitatively relevant also for the beam 2, which does the counter clockwise 
rotation. 
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Figure 1. The LHC magnet structure close to the ATLAS Interaction Point. 



The ALFA experiment consists of four detector stations placed symmet- 
rically with respect to the ATLAS IP at 237.4 m and 241.5 m. In each 
station there are two roman pot devices, which allow to insert the position 
sensitive and triggering detectors vertically into the beam pipe. Two sta- 
tions are needed at each side to be able to measure not only the scattered 
proton trajectory position, but also its direction (elevation angles). 

An important point is to understand the dependence of the scattered 
proton trajectory position in the detector on its energy and momentum. 
It is illustrated in Figs [2a| and [2b] for early and nominal high /3* optics, 
respectively. For both optics settings the impact of the pF'-momentum 
component at the IP on the proton position in the detector station is much 
greater than that due to p^ . The higher is the proton energy loss, AE, 
the higher is its deflection towards the machine centre. However, due to the 
differences in the LHC optics for both tunes this deflection is larger in the 
case of the nominal high (3* . 

Naturally, not all scattered protons can be measured in the ALFA de- 
tectors. Such proton can be too close to the beam to be detected or it can 
hit the LHC elements (a collimator, the beam pipe) in front of the ALFA 
station. The geometric acceptance, shown in Fig. [3] for both optics settings, 
is defined as a ratio of the number of protons of a given energy loss (A£') 
and transverse momentum (pt) that crossed the active detector area to the 
total number of the scattered protons having E and px- In the calculations 
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Figure 2. The positions of protons with different energy loss (Ai?) and transverse 
momentum (pt) at the first ALFA station for the two different LHC optics settings. 
The solid lines mark the beam pipe aperture and the ALFA detector active area. 



the following factors were taken into account: the beam properties at the 
IP, the beam chamber and the detector geometries, the distance between 
the detector edge and the beam centre. This distance was set 4 and 1.5 mm 
- the values expected for the early and nominal high /3* runs. Values of 
other parameters are listed in Tabled) 

As can be observed, the region of acceptance above 80% is limited by 
AE < 0.6 TeV and 160 MeV/c < pr < 200 MeV/c for early high /3*. If the 
requirement on the acceptance value is lowered to 60% then the range of the 
accepted transverse momentum values gets slightly larger. In the case of the 
nominal high /3* the high acceptance region is a triangle-like and spans a 
bit larger range of the proton transverse momentum. Moreover, it is worth 
mentioning that high /3* runs will have very low both instantaneous and 
integrated luminosities. This implies that only a few events containing par- 
ticles particles with high energy loss are expected to be observed. Therefore, 
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Figures. The geometrical acceptance of the ALFA detector as a function of the 
proton energy loss (A£') and its transverse momentum (pr) for two LHC settings. 
The distance between the beam centre and the detector edge was set to 4 mm for 
early high /?* and for 1.5 mm for nominal high /3*. 



Table 1. The LHC beam parameters at the ATLAS IP for early and nominal 
high (3* runs. 



Parameter Unit Early High /3* Nominal High /3* 






'pt 



mm 
/irad 

MeV/c 



0.3 
3.33 

n.7 



0.612 

0.233 

1.6 



the most important factor is the accepted pT range for small /S.E. 

The minimum transverse momentum of protons that can be registered 
depends on the distance between the detector edge and the beam centre. 
This is demonstrated in Fig. H] where the scattered proton transverse mo- 
mentum spectra are shown for different distances and both machine settings. 
Clearly, the smaller the distance between the beam and the detector, the 
smaller is the limiting value of the accepted proton's px- This is particu- 
larly important for the elastic scattering measurement where the possibility 
of reaching as small pT values as possible is crucial. 
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Figure 4. The influence of the distance between the detector edge and the beam 
centre on the accepted proton transverse momentum. 



3. Transport Parametrisation 

Anticipated positions of protons at the detector for a given momentum 
can be obtained from the transport simulation. One can also prepare a 
look-up table and interpolate the positions. The disadvantage of the above 
methods is either a long calculation time as in the first case or very extensive 
use of the computer storage (as size of the looking table grows with the 
required precision) in the second. The alternative is to use the transport 
parameterisation, which is very fast and requires only very little memory 
space. Moreover, the parameterisation provides an analytical representation 
of the proton position and momentum. This idea was first proposed in [7]. 

The LHC magnetic structure in the vicinity of the ATLAS detector is 
described only by the drift spaces, the dipole and the quadrupole magnets 
(c/. Fig. [1] or [8J). Therefore, a linear transport approximation can be 
applied to describe the scattered proton transport. 

For a transverse variable C € {j;, y, j;', y'} the transport can be effectively 
described by the following equations: 



Ac^^Bc^- xip + Cc_ ■ yip + -De • ^ip + ^^ ' ^ip 



(1) 



where A,^, . . . ,H(^ are the polynomials in the reduced energy loss (^ 



LHC'High'Beta'Runs'Transport'and'UnJolding'Methods printed on July 12, 20111 



/\E/Ebeam) of rank A;^ . . . , /c^.: 






kHr 



^c = E«c,n-r, ■■■, H^ = J2k,n-c. 



n=0 



n=0 



The absence of magnets with multipole field expansion moments higher 
than the quadrupole one imphes that the horizontal trajectory position (di- 
rection) does not depend on the vertical momentum component nor the 
vertical vertex coordinate, and vice versa. The best description of the scat- 
tered proton transport for the ALFA detectors is given by: 



t^Ax ^Bx ^Ex 

E «x,n • ^" + E bx,n • C" • a;iP + E ex,n ' C" " x[p, 
n=0 n=0 ra=0 

kCy kpy 



y 



E Cy^n • r • 2/IP + E fy,n ■ T " 2/ff, 
n=0 ra=0 

E ax',n • C" + E K',n • C" • a;ip + E ex',n • C" • a^ip, 

n=0 ra=0 71=0 

E S',„-r-?/ip+ E /j/'.n-r-yip, 

n=0 n=0 



(2) 



The exact values of the coefficients were found by fitting Eq. ([T]) to the 
MAD-X PTC results. This parameterisation was validated with an inde- 
pendent single diffractive event sample generated with Pythia 6.4 Monte 
Carlo [9]. In the simulation the interaction vertex position was smeared 
appropriately for the discussed LHC tunes (see Tab. [1]). The parametrisa- 
tion uncertainty was evaluated comparing the results obtained with Eq. ([2]) 
to those of the MAD-X PTC. Results of this comparison are presented in 
Table [2j The differences between the MAD-X PTC and the parametrisation 
results are much smaller than the detector resolution (30 fim). This con- 
firms that the parametrisation provides a good description of the scattered 
proton trajectory positions and elevation angles at the detectors positions. 



4. Unfolding Procedure 

The procedure of inferring the scattered proton momentum on the basis 
of the detector measurements is called unfolding. This can be done in 
various ways, here is is performed by means of the minimisation of the 
following x^ function: 



X\p) 



(xf 



2;UP, 



at. 



-+ 



y? 



yi(p) 



-+ 



X. 



D 



2;2(p)) (; 



y^ 



erf, 



cr:t 



2/2 IP 



cff, 



SLHC'High'Beta'Runs'Transport'and'Unfolding'Methods printed on July 12, 2011 



Table 2. The parameterisation method uncertainty measured as the RMS of the 
difference between values given by MAD-X PTC and parameterisation equations 
for the early and nominal high /3*. 



Variable Unit Early High /?* Nominal High /3* 



Ax 


nm 


45.1 


37.9 


Ay 


nm 


50.4 


22.9 


Ax' 


nrad 


35.7 


4.7 


Ay' 


nrad 


9.1 


10.3 



where (xf,yf) denote the coordinates of the scattered proton trajectory 
position measured by the first station and (xi(p),yi(p)) - the coordinates 
calculated using the transport parametrisation for a proton with momentum 
p. The variables X2, yf', 3^2 (p) and 1/2 (p) refer to the positions at the second 
station. The parametrisation does not describe correctly neither the losses 
of particles due to the beam pipe nor the collimators apertures, but this 
drawback is of a minor importance for solving the unfolding problem. 

Results of the unfolding procedure are presented in Fig. [5l The corre- 
lations between the reconstructed and the generated value of the reduced 
energy loss and the transverse momentum components of the scattered pro- 
ton are shown for the early high /3* runs. For comparison the horizontal 
transverse momentum unfolding results are added. This figure presents also 
the influence of such experimental factors as the vertex smearing and the de- 
tector spatial resolution [ax = ay = 30 fiTii is assumed). The leftmost plots 
in Fig. [5] confirm the quality of the parametrisation and the correctness of 
the unfolding procedure. If particular experimental effects are taken into 
account the unfolding accuracy worsens. Due to the chosen optics mode, 
the momentum reconstruction works better for the y coordinate. The re- 
sults show that in the early high f3* case the vertex smearing has a major 
impact on the horizontal momentum reconstruction error whereas in the 
nominal high /3* the leading effect is governed by the detector resolution. 
The difference in the horizontal momentum resolutions are due to the x- 
vertex spread. This is shown in Fig. [6l where a correlation between the px 
unfolding error and the horizontal vertex coordinate is plotted. It clearly 
shows that the x-vertex spread has a much greater effect for the nominal 
high /?* reconstruction accuracy. 
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Figure 5. Correlations between the true and unfolded values of reduced energy 
loss and transverse momentum components. The left column shows the ideal case 
where the measurement is perfect, whereas the influence of experimental effects, 
the vertex smearing and the detector resolution, are presented in the middle and 
right column, respectively. 
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Figure 6. Left: the correlation between the horizontal momentum reconstruction 
error and the horizontal vertex coordinate. Right: the reduced energy loss recon- 
struction resolution. 



It is interesting to study what precision can be obtained for measure- 
ments of the diffractively scattered protons. This would yield information on 
possible measurements of diffraction that can be performed with the ALFA 
detectors. The single diffractive dissociation events were generated with 
Pythia. The final state forward protons were transported to the ALFA 
stations and the actual positions of the trajectories were smeared according 
to the detector spatial resolution. 

Then, the unfolding procedure was performed and its errors were esti- 
mated by fitting the distribution of the difference between the generated 
value and the one obtained from the unfolding with the Gaussian distribu- 
tion. These fits were performed for different bins of ^ (Fig. [6] right), px 
and Py (Fig. [7]). For comparison the momentum resolutions were calculated 
for the elastically scattered protons. In this case, the unfolding procedure 
was slightly different, because the x^ minimisation was performed only in 
Px and Py variables as the proton energy is fixed for such events. 

One immediately notices that the early high /3* offers better opportu- 
nities for obtaining diffractively scattered proton kinematics. In that case 
the energy loss reconstruction error is about 40 GeV which compares to 
160 GeV for the nominal high /3*. One should notice large differences be- 
tween the resolutions in the horizontal and the vertical directions. This is a 
direct consequence of the parallel-to-point focusing optics feature dedicated 
to a precise measurement of small scattering angles in the vertical direction. 

In the case of elastic scattering the situation is opposite: the horizontal 
momentum reconstruction for the nominal high /3* optics has much better 
resolution than that for the early high /3* one. The vertical momentum is 
reconstructed very accurately and its reconstruction resolution is about 0.3 
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Figure 7. The horizontal (left) and vertical (right) momentum resolution for elas- 
tically and diffractively scattered protons for early and nominal high /3* optics. 

MeV/c for both cases. This is because the nominal high (3* optics is specially 
designed for very precise elastic scattering measurements. However, one has 
to remember that in the case of all high /3* runs, the beam angular spread 
can be a serious limiting factor (c/. Tab. [TJ. For example, for the early 
high /3* tune the beam transverse momentum spread is about 20 times 
bigger than the py reconstruction resolution. 



5. Summary 

The transport of the elastically and diffractively scattered protons through 
the LHC magnetic lattice for the high /3* optics settings was described. The 
geometrical acceptance for the two most probable ALFA run settings: the 
early high /3* (E'beam = 3.5 TeV, /3* = 90 m, e* = 2.5 /^m-rad) and the 
nominal high f3* (ii'beam = 7 TeV, /3* = 2625 m, e* = 1 /im-rad) was pre- 
sented. Studies of the transverse momentum acceptance as a function of 
the detector edge distance from the beam centre reveals that the discussed 
/3* settings are the most desirable for the luminosity determination because 
the region of much lower pT values can be accessed. 

The transport parameterisation was introduced as the fastest simulation 
method of the particle transport through the LHC magnetic lattice, which 
delivers an analytical representation of the scattered proton position at the 
detector stations. The accuracy of this method was shown to be much better 
than the assumed spatial resolution of the detector. 

Finally, the unfolding method was presented as a procedure to extract 
the scattered proton energy and momentum at the Interaction Point from 
its trajectory measurement at the forward detectors. The example of the 
ALFA case shows that the proton energy can be reconstructed with preci- 
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sion of about 40 GeV in the case of early high /3* and about 160 GeV for 
the second discussed setting. The scattered proton momentum reconstruc- 
tion precision is dominated by its horizontal component resolution and is 
about 0.15 GeV/c for both optics. In the elastic scattering case the momen- 
tum resolution is still dominated by its horizontal component reconstruction 
resolution, but the nominal high (3* setting allows for twenty times better 
precision than the early high /3* one. 
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